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^ ■ Abstract. Relationship between quantum shell structure and classical periodic orbits is briefly reviewed 

■^r ' on the basis of semi-classical trace formula. Using the spheroidal cavity model, it is shown that three- 

I dimensional periodic orbits, which are born out of bifurcation of planar orbits at large prolate deformations, 

^— ^ , generate the superdeformed shell structure. 



Introduction 



> 

OO ' Existence of superdeformed (SD) nuclei is often explained in terms of the SD magic numbers for the harmonic- 

. oscillator (HO) potential with axis ratio 2:1. It appears, however, that we need a more general explanation not 
' restricted to the HO potential, since, up to now, more than 200 SD bands have been found in various regions 
of nuclear chart and their shapes in general deviates from the 2:1 shape to some extent. In this talk, we shall 
discuss the mechanism how and the reason why the SD shell structure emerges. The major tool for this purpose 
is the trace formula, which is the central formula in the semiclassical periodic-orbit (PO) theory and provides 
a link between quantum shell structure and classical periodic orbits in the mean field. Here, shell structure is 
^ I defined as regular oscillation in the single-particle level density coarse-grained to a certain energy resolution. 
P5 . An example of coarse- graining for the well-known axially symmetric HO model is displayed in Fig. 1. 

In this talk, we discuss the spheroidal cavity model, since, in contrast to the HO model, this model is very 
rich in periodic orbits; it is an ideal model for exhibiting the presence of various kinds of periodic orbit and 
' their bifurcations. We present both Fourier transforms of quantum spectra and semiclassical calculations based 
H \ on the PO theory, and identify classical periodic orbits responsible for emergence of the SD shell structure. 

The result clearly shows that three-dimensional (3D) periodic orbits, that are absent in spherical and normal 
deformed systems and are born out of bifurcations of planar orbits, generate a new shell structure at large 
prolate deformations, which may be called "the SD shell structure." They continue to exist for a wide range 
of deformation, once they are born. 

The PO theory provides a basic tool to get a deeper understanding of microscopic origin of symmetry 
breaking in the mean field. It sheds light, in addition to the stability of the SD nuclei, on the reason of prolate 
dominance in normal deformed nuclei, on the origin of left-right asymmetric shapes, etc. It is useful for finite 
many-Fermion systems covering such different areas as nuclei, metallic clusters, quantum dots, etc. In this 
talk, we shall also touch upon such applications of the PO theory. 

Level Bunching and Trace Formula 

For the axially symmetric HO potential, the following two conditions coincide: 
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FIGURE 1. Oscillating level density as a function of axis ratio rj for the axially symmetric HO potential. Bright and 
dark regions correspond to high and low level densities, respectively. The 3D periodic orbits exist only at the 2:1 shape. 



1) Axis ratio {x"^) : {z^) = 1:2 

2) Frequency ratio lDx '■ LiJz = "2. '■ 1 , 

but they are different in general. Condition 2) is nothing but the PO condition, and possesses a more general 
significance than condition 1). We can examine this point as follows. For any integrable Hamiltonian system, 
we can introduce action and angle variables (J, 6) which satisfy the canonical equations of motion, 

ri M 

^ = 57 - ^(i)^ (1) 

and the energy E can be quantized by the EBK (Einstein-Brillouin-KcUer) quantization condition: 

En = H(I) with I = h{n+\cx), (2) 

where n represents a set of quantum numbers, n = {n\,n2, na) with rii = 0, 1, 2, • • •, and a the Maslov indices. 
We see that level degeneracy occurs when 

K+Ar^ -Er,= H {I + AJ) - H{I) 

= hu) ■ An 

= fkviAm + fiu!2An2 + hwsAn^ 

= 0, (3) 



i.e., when uji : uj2 '■ W3 are in rational ratios. This is just the condition for the classical orbit to be periodic, 
and discussed in detail in the textbook of Bohr and Mottelson [1]. 

On the basis of the semiclassical PO theory, we can examine, in a more general way, the decisive role of 
periodic orbits as origin of level bunching. According to this theory (see, e.g., [2] for a review), the level density 
g{E) is given by a sum of the average part g{E) and the oscillation part Sg{E) as 

g{E) = Y,S{E-E^) 

n 

c^g{E)+5g{E) 

= ff(£;)+5^A„cosQ5„(E)-|/Xa) , (4) 

where Sa{E) denotes the action of the periodic orbit a, and is a phase related with the Maslov index. This 
equation is called "trace formula" and provides a link between quantum shell structure and classical periodic 
orbits in the mean field. There is a complementarity between the energy resolution and periods of classical 
orbits, so that, for the purpose of understanding the origin of regular oscillation patterns in the smoothed 
single-particle level density (shell structure), we need only short orbits in the sum over a. 



Periodic Orbits and Shell Structure in the Cavity Model 

Let us consider the cavity model, which may be regarded as a simplified model of Woods-Saxon potential 

for heavy nuclei. In fact their basic patterns of shell structure are similar with each other. Certainly, the 
spin-orbit term shifts the magic numbers, but it does not destroy the valley-ridge structure discussed below. 
One can confirm these points by comparing Figs. 2, 3 and 4. 
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FIGURE 2. Oscillating level density as a function of 

cixis ratio rj for the WS potential without the spin-orbit 
term. Bright and dark regions correspond to high and 
low level densities, respectively. The calculation was done 
following the procedure described on p. 593 of Ref. [1]. 




FIGURE 3. Oscillating level density as a function of 

axis ratio 77 for the WS potential. Bright and dark regions 
correspond to high and low level densities, respectively. 
Here, the spin-orbit term with vu = —0.12 is added to 
the Hamiltonian used in Fig. 2. 
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FIGURE 4. Ridge-valley structure in the oscillating level density for the spheroidal cavity model. Bright and dark 

regions correspond to high (ridge) and low (valley) values of the oscillating level density, respectively. Constant action 
lines running along valleys are shown for typical periodic orbits by dashed, thin-solid, and thick-solid curves. Note 
that they arc indicated as representatives of the meridian-plane orbits, the equatorial-plane orbits, and the 3D orbits, 
respectively: There are several families of periodic orbits with similar lengths, and constant-action lines of them also 
behave in the same way as those shown here, see Refs. [4] and [10] for details. 



Role of periodic orbits for shell structure in this model was originally studied by Balian and Bloch [3] , and 
has been discussed in the investigations of 

1) the reason of prolate dominance in nuclear shape [4-6] 

2) the supershell effects in metallic clusters [7] 

3) the origin of left-right asymmetric shapes [8,9], etc. 

For cavity models, the energy and the momentum are simply related as 

E = p^/2m, p = hk, (5) 
and the action integral is proportional to the length La, 

Sa = f p - dx = hkLa- (6) 

Accordingly, the trace formula for the level density can be written as 

9{k) gik) + ^AaCos(ji:La-^IJ,a^ ■ (7) 

a 

It can be easily confirmed that only short orbits contribute to the level density coarse-grained in energy. Let 
us Fourier transform the level density 



F{L) = j dke'''^g{k) 

^Y.A^5{L-L^). (8) 

a 

This equation indicates that peaks wiU show up at lengths of periodic orbits a, which may be called "length 
spectrum" . Now, the orbit lengths change when the deformation parameter rj varies. Let us then consider the 
oscillating level density as a function of r}, 

6g{k, r?) ^ ^ Aa{k, rj) cos (^kLa{v) - -^Ma) • (9) 

a 

Prom this formula, we see that, if a few orbits dominate in the sum, the valley-ridge structure on the {k, rj) 
plane will be determined by the constant action lines, 

kLa{rj) = const, (10) 

of these dominant orbits [4,10]. In fact, we see in Fig. 4 that the valley-ridge structure is well explained in 
terms of three kinds of short periodic orbit. 

Bifurcations 

As illustrated in Figs. 5 and 6, when the axis ratio r] of the spheroidal cavity reaches \/2, the butterfly-shaped 
planar orbits emerge on the meridian plane through bifurcations of the linear orbit along the short diameter. 




FIGURE 5. Birth of a butterfly shaped orbit from the short diameter through bifurcation at rj = \/2 on the meridian 
plane. This figure illustrates a representative orbit among a continuous family of orbits possessing the same length. 




FIGURE 6. Birth of a 3D orbit from the star-shaped orbit on the equatorial plane through bifurcation at r? ~ 1.62. 
This figure illustrates a representative orbit among a continuous family of orbits possessing the same length. 
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FIGURE 7. Left-hand part: Absolute values of the Fourier amplitudes, \F{L,r])\, of the level density, drawn as a map 

on the (L, Tj) plane. Brightness is proportional to the magnitude. Right-hand part: Lengths of various classical orbits as 
functions of axis ratio rj. The 3D orbits responsible for the SD shell structure are denoted as (4,2,1), (5,2,1), (6,2,1), etc., 
and their lengths are plotted by thick-solid curves, while those of the equatorial-plane orbits and of the meridian-plane 
orbits are plotted by thin-solid and dashed curves, respectively. ©The bifurcation points of the 3D orbits of interest are 
indicated by filled circles. (Filled triangles indicate bifurcation points of 3D orbits responsible for hyperdeformed shell 
structure not discussed here.) See Refs. [10] and [13] for details. 



When 1] further increases, 3D orbits emerge at rj ~ 1.62 through bifurcation of the five-point star shaped orbits 
on the equatorial plane. Likewise, other 3D orbits appear a.t r] = -\/3 through bifurcation of second repetitions 
of the triangular orbits on the equatorial plane, . . . , etc. Note that the figure illustrates representative orbits 
only. In fact, each a in the trace formula (9) represents a continuous family of orbits with the same topology 
possessing the same values of action (length). In contrast to the HO potential, these 3D orbits continue to 
exist, once they appear through the bifurcations. 

Peaks in the Fourier transform (8) of the level density will follow the variations of orbit lengths with rj. 
Thus, we can draw a map of the Fourier amplitudes on the {L, rj) plane, 

F(L,r,) = ^i„<5(i-L„(77)). (11) 

a 

In Fig. 7, the Fourier amplitudes are compared with lengths of classical periodic orbits. This figure exhibits 
a beautiful quantum-classical correspondence. Furthermore, by comparing the bright regions in the left-side 
figure with the bifurcation points indicated in the right-side figure, we find significant enhancement of the shell 
structure amplitudes just on the right-hand side of the bifurcation points. 

Unfortunately, the amplitude Aa{k^ri) in the trace formula (9) diverges at the critical point of deformation 
T] where the orbit bifurcation takes place. This is because the stationary phase approximation used in the 
standard semiclassical PO theory breaks down there. Thus, the standard trace formula is unable to describe 
the enhancement phenomena seen in Fig. 7. To overcome this difficulty, in recent years, we have developed a new 
semiclassical approximation scheme, called an improved stationary phase approximation, and derived a new 
trace formula free from such divergence [11-13]. A numerical example obtained by this approach is presented 
in Fig. 8. We see that the basic pattern of oscillation in the quantum level density at large deformation is 
nicely reproduced by the semiclassical calculation using the new trace formula. In this way, we have confirmed 
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FIGURE 8. Comparison between the oscillating level densities evaluated by the semiclassical trace formula (dotted 
lines denoted ISPM) and those obtained by quantum mechanics (solid lines denoted Quantum). They are plotted for 
several values of r; as function of kR, where R denotes the radius of the cavity in the spherical limit (the volume of 
cavity is conserved during variation of r]) . 



[10,12,13] that, in the region of large prolate deformation with axis ratio t] > 1.62 (which corresponds to the 
ordinary deformation parameter S > 0.44), the major pattern of the oscillating level density is determined by 
contributions from the bifurcated 3D orbits. 



Conclusion 



The 3D periodic orbits generate a new shell structure at large prolate deformations. Wc may call this shell 
structure "SD shell structure." These 3D orbits are born out of bifurcations of planar orbits in the equatorial 
plane, and they play dominant roles immediately after the bifurcations. Thus, the SD shell structure is 
a beautiful example of emergence of new structure through bifurcation, and may be regarded as quantum 
manifestation of classical bifurcation. 
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